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Recent experiments have unveiled important properties of the ground state of the elusive heavy 
fermion URu 2 Si 2 . While tetragonal symmetry-breaking was reported below the hidden-order (HO) 
transition at Tho ~ 17-5 K, time-reversal symmetry-breaking was observed below the superconduct¬ 
ing transition temperature Tc < Tho- Although the latter results have been used to argue in favor 
of a chiral d + id superconducting state, such an order parameter is incompatible with broken tetrag¬ 
onal symmetry. Here, we employ a phenomenological model to investigate the properties of a chiral 
superconducting state that develops inside the hidden-order phase. In this case, there are actually 
two superconducting transition temperatures: while marks a normal-state to superconducting 
transition, Tf < Tc signals a superconducting-to-superconducting transition in which time-reversal 
symmetry is broken. In the phase T* < T < Tc, the low-energy density of states p{to) is enhanced 
due to the crossing of two nodal lines, giving rise to an unusual a;log(a;) dependence of p{u)), which 
is manifested in several thermodynamic properties. We also investigate the emergence of a soft 
amplitude gap mode near T*. In contrast to the usual amplitude mode near a regular normal-state 
to superconducting transition, this mode becomes soft near a superconducting-to-superconducting 
transition, which in principle allows for its detection by Raman spectroscopy. Finally, we investigate 
the impact of twin domains on the anisotropic properties of the superconducting state, and propose 
experiments in mechanically strained samples to explore the interplay between hidden order and 
superconductivity in URu 2 Si 2 . 


I. INTRODUCTION 

The nature of the ground state of the body-centered 
tetragonal compound URu2Si2 remains one of the most 
elusive problems in heavy fermion physics. At the 
HO transition temperature Tho ~ 17.5 K, the sys¬ 
tem displays a sharp specific heat anomaly characteris¬ 
tic of a second-order phase transitioi>i"— . However, de¬ 
spite nearly thirty years of intense research, the bro¬ 
ken symmetries of this hidden-order phase remain a 
widely debated topic. On the theory front, many 
different order parameters have been proposed to ex¬ 
plain the hidden-order state, from high-rank multiple 
ordersi^— to exotic states involving the hybridization 
of localized and itinerant state&ii^— . On the experi¬ 
mental front, recent measurements have provided im¬ 
portant pieces for this long-standing puzzle. For in¬ 
stance, torque magnetometry^ii^^, x-ray difFractioi>^, 
and elasto-resistance measurements^ have reported a 
tetragonal {C4) to orthorhombic {C2) transition simulta¬ 
neously to Tho , manifested by the inequivalence between 
the [110] and [l —1 O] directions of the crystallo¬ 
graphic unit cell (see Fig. |T]). Whether this is the only 
broken symmetry in the hidden-order state - in which 
case it would be classified as a nematic state - remains 
an open issu o^^i^^ . For instance, quantum oscillation 
measurements also suggest a translational symmetry¬ 
breaking along the c-axis^, and neutron elastic scatter¬ 
ing experiments favor a rank 5 multipole order in the 
system^. 

In comparison to the hidden-order phase, the super¬ 
conducting state of URu2Si2 that appears at Tc ~ 1.5K 
has received rather less attention. Understanding its na¬ 
ture is relevant not only within the bigger picture of su¬ 


perconductivity in heavy-fermion compounds, but also 
as a potential tool to probe the properties of the HO 
state^i^, since superconductivity develops well below 
Tho and disappears when the HO phase is suppressed by 
external pressure^i^. An interesting proposal based on 
recent angle-resolved specific heat^ and thermal conduc¬ 
tivity data^, which reported indirect evidence for point 
and line nodes, is of a chiral d-wave superconducting (SC) 
state described by the order parameter—: 

A(fc) = Aosmy ( sm-y-^±^sm-. (1) 

Such a SC state breaks time-reversal symmetry, which 
seems to be in agreement with magnetic susceptibility^ 
and recent Kerr effect measurements^. However, this 
order parameter manifestly preserves C4 symmetry, 
whereas the HO state, from which SC develops, breaks 
tetragonal symmetry. 

In this paper, we use a phenomenological model to 
reconcile the proposal of a chiral d-wave state with the 
experimental observations of a HO state that breaks C4 
symmetry. While a phenomenological approach leaves 
aside the issue of the microscopic mechanisms involved 
in the formation of these phases, it allows for general 
conclusions to be drawn regardless of one’s favorite or¬ 
der parameter for the HO phase - as long as it accounts 
for C4 symmetry breaking. As a result, it provides gen¬ 
eral benchmarks that must be satisfied if indeed the HO 
state breaks C4 symmetry and the SC state is chiral. One 
obvious consequence from the fact that SC develops in a 
C2 background, as pointed out in different scenarios^^— , 
is that the chiral superconducting transition is actually 
split into two superconducting transitions Tc and T*. 
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While at Tc the system first becomes superconducting, 
at T* < Tc time-reversal symmetry (TRS) is broken, sig¬ 
naling a SC-SC transition. Our focus in this paper is on 
the thermodynamic properties of these two distinct su¬ 
perconducting phases. We find that, while in the regime 
T < T* the nodal quasi-particle density of states p(w) 
depends linearly on w for low-energies, reflecting the pres¬ 
ence of point and line nodes, in the regime T* < T < Tc 
the density of states acquires an unusual log-dependence 
p{ui) ^ wlogw due to the crossing between two nodal 
lines. Such a behavior leaves signatures in several ther¬ 
modynamic quantities, such as the specific heat and the 
penetration depth. 

Furthermore, we investigate how the anisotropies of 
the SC state in both regimes - namely T < T* and 
T* < T < Tc - are affected by the presence of twin 
domains with different C2 orientational order. By cal¬ 
culating the angular dependence of the specific heat in 
the presence of a magnetic field, we find that a twinned 
sample and a sample with no C4 broken symmetry would 
display nearly identical behaviors. This helps to reconcile 
some of the experimental results that led to the proposal 
of a chiral SC state, as in Eq. o, with the experimen¬ 
tal results that found tetragonal symmetry-breaking in 
the HO phase. We propose experiments in mechanically 
detwinned samples to unveil the intrinsic anisotropies of 
the SC state. 

We also find an unusual behavior for the collective SC 
modes near T*. In a common normal-state to SC transi¬ 
tion, the soft amplitude SC gap-mode falls into the con¬ 
tinuum and is strongly damped. However, because T* 
is a SC-SC transition, the corresponding soft amplitude 
mode develops in the background of a superconducting 
quasi-particle spectrum. Despite the presence of nodal 
quasi-particles in this spectrum, which promote under¬ 
damping of the soft amplitude mode, we argue that this 
collective mode may still be observed in the excitation 
spectrum, as measured by Raman scattering. 

The paper is organized as follows: in Section H we 
present our phenomenological model for the SC state in¬ 
side the HO phase and discuss its nodal quasi-particle 
spectrum for different temperature regimes. In Sec¬ 
tion HI we study the impact of twin domains on the 
anisotropic thermodynamic properties of the SC state, 
with particular emphasis on the angle-dependent specific 
heat. Section IV is devoted to the investigation of the 
collective modes of the SC state. Concluding remarks 
follow in Section V. 


II. PHENOMENOLOGICAL MODEL 

A. Superconducting free energy 

The SC order parameter in Eq. m belongs to the 
Eg irreducible representation of the tetragonal point 
group. Its degeneracy stems from the two-dimensionality 
of the Eg representation. However, because the HO state 


breaks the C4 symmetry, the SC order parameter must 
be modified to reflect the new orthorhombic point group 
symmetry of the system. Using as a starting point Eq. 
m, we can describe the chiral d-wave SC order parameter 
inside the orthorhombic HO phase as a two-component 
order parameter: 


A{k) = Ai sin — sm —- 

, iO A . ky . , 

+e Aasmysm -^ , (2) 

where 0 is the phase difference between the two com¬ 
ponents of the order parameter, and are real order 
parameters. TRS is kept intact only if 0 = 0,7r. The 
absolute value of the gap function is then: 


kj; -j- ky 


ky- kj, 


I A |2 • 2 A 2 * 2 “T Ay A 2 • 2 

I A| = sm — A^ sm - - + Ai, sm - - 


^ A A Ai . kx ky 

-|-2Ai A2 cos 0 sin-- —- sm 


■2 

kx — k. 


( 3 ) 


Thus, tetragonal symmetry requires 0 = and 

Ai = A2, but in the HO phase this is not the case. 

To proceed, we need to establish how the tetragonal 
symmetry is broken in the hidden order phase. There 
are two possibilities: either the horizontal [l 0 O] and 
vertical [O 1 O] directions become inequivalent, in which 
case the C2 order parameter belongs to the Big irre¬ 
ducible representation, or the diagonal [l 1 O] and anti¬ 
diagonal [1 —1 0] directions become inequivalent, in 
which case the C2 order parameter belongs to the B^g 
irreducible representation. Torque magnetometry^^, x- 
ray diffraction^, and elasto-resistance measurements^! 
all agree on the second scenario, which is schematically 
shown in Fig. [T] 
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Figure 1. The basal plane across the HO phase transition. 
Left: above the HO phase, the crystal structure (here centered 
at the U atoms) is body centered tetragonal. Right: The HO 
phase breaks the fourfold rotational symmetry on the ah plane 
by making the [l 1 O] and [l —1 O] directions (parallel to 
the Ru nearest neighbor directions) inequivalent. The crystal 
structure becomes base centered orthorhombic. 


Thus, we can now write down the phenomenological 
Ginzburg-Landau model for the superconducting degrees 
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of freedom inside the HO phased: 

Fsc = ^(A? + A^) + |(At + A4) 

+ ^A2a 2 (/3 + acos26») - I (A? - A^) , (4) 


A? = - 


F — — 

m 1 n — 


a — Tj 

U 


{a - v f 

4 m 


A2 = 0 


( 6 ) 


Here, a = qq (T — T^fi) and the order parameter ^ 0 
describes the symmetry-breaking inside the HO state. 
We refrain from calling it the HO order parameter, since 
it is unclear whether other symmetries are also broken 
in the HO phase. In any case, because Tho ^ 7 c> we 
consider ry to be constant on the temperature regime in 
which an expansion in powers of the superconducting or¬ 
der parameter is allowed. Without loss of generality, we 
assume in this section that ry > 0 - in the next section, 
where we consider the effects of domains, this assumption 
is no longer valid. 

For the SC state to be a chiral d-wave, a must be 
positive; furthermore, in order for both order parame¬ 
ter components to coexist, and for the free energy to be 
bounded, we must also have u > |/3 — a\. With this con¬ 
straints, minimization of the free energy with respect to 
6 gives always 9 = ±7r/2, as long as both Ai and A2 are 
simultaneously non-zero. 


A/ 



Figure 2 . The SC order parameters in different phases. Ai 
becomes nonzero when T < Tc (SCi phase), and A2 becomes 
nonzero at T < T* (SC2 phase). The system, therefore, has 
two different SC phases. The phase difference between Ai 
and A2 is tt/2, i.e. time reversal symmetry is broken only in 
the SC2 phase. 


Minimization of the free energy with respect to Ai and 
A2 reveal two distinct regimes as shown in Fig [2 defined 
by the two transition temperatures: 


Tc — Tcfl + 


Oo 


T: = T, - 


2ur] 


ao{u — P + a) 


( 5 ) 


Here, we assume that a = ao(T — when the tem¬ 
perature is close to Tcfi. For T* < T < (SCi phase), 
the free energy is minimized by enforcing only one of the 
SC components to be non-zero: 


On the other hand, for T < T* (SC2 phase), the mini¬ 
mum of the free energy corresponds to the condensation 
of the two components with a 7 r /2 phase difference: 


9 = 'K 12 , 



u-\-l 3 — a~^u — j 3 -\-a ’ 
a rj 

u + P — a u — P + a ^ 


2(u + p — a) 2{u — P + a) 


( 7 ) 


Because the free energy changes smoothly, both tran¬ 
sitions are second-order. At Tc, the system becomes a 
single-component SC, whereas at T*, TRS is broken, and 
the system becomes a two-component SC. The anisotropy 
induced in Ai and A2 is proportional to ry, as expected, 
since this order parameter manifestly breaks the tetrag¬ 
onal symmetry. The splitting between Tc and T* is also 
proportional to | 7 y|. 

It is important to discuss the relationship between the 
orthorhombic distortion S in the C2 HO phase and the 
electronic anisotropy order parameter rj introduced here. 
Symmetry arguments enforce them to be proportional 
to each other, i.e. {S) oc (?y). However, the fact that S 
is small in the HO phase {6 ~ 10 “® as measured by x- 
ray diffraction^) does not imply that rj is also necessarily 
small. If indeed the anisotropy is electronically driven, as 
suggested by the elasto-resistance measurements^!, one 
would expect jy to be sizable even if 6 is small. This is the 
case, for instance, in optimally-doped iron-based super¬ 
conductors, where the in-plane resistivity anisotropy is 
Pb/Pa ^ 1-5 even though 6 ~ 10 “'* (for reviews, se o^'*'^^ ). 


B. Nodal quasi-particle density of states 

Having established the existence of two SC states in 
URu2Si2, we now discuss their thermodynamic prop¬ 
erties. At low temperatures, they are determined by 
the low-energy properties of the quasi-particle density 
of states (DOS) p{u}). For instance, the specific heat C, 
the penetration depth Ap, and the spin-lattice relaxation 
rate l(TiT are given by: 

^ oc y dujpiuj) {Pujf 

V" - V.o J dujpiuj) 

(TiT)"^ oc j dup^ (w) 


( 8 ) 
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At low energies, the DOS is determined by the dis¬ 
persion of the nodal quasi-particles. The latter are re¬ 
markably different for the two SC states. For the higher- 
temperature SC state at T* <T < Tc, which we denote 
SCi, only one of the gap components is condensed and 
the gap function contains horizontal zeros at kz = 0, ± 27 r 
and vertical zeros at kx = —ky (if Ai ^ 0) or at 
kx = ky (if A2 ^ 0 ). On the other hand, for the lower- 
temperature SC state at T < T*, which we denote SC2, 
the gap vanishes along the plane kz = 0 ,± 27 r and along 
the line kx = ky = 0. 


To proceed, we need to establish whether these gap ze¬ 
ros cross the Fermi surface, forming nodal quasi-particles. 
The Fermi surface of URu2Si2 is remarkably complex: 
first-principle calculations in the paramagnetic phase 
show the existence of hole-like ellipsoids centered at F = 
( 0 , 0 , 0 ) and Z = ( 0 , 0 , 27 r), and electron-like ellipsoids 
centered at M = (tt, tt, 0 ) 27 ^^ Although it is unclear 
how the hybridization with the 5 / U states affects this 
Fermi surface configuration, or whether there is an ad¬ 
ditional folding of the band structure along the (0,0, 27 r) 
momentum, as in the magnetically ordered phase, quan¬ 
tum oscillation measurement o"^°'^^~ — seem to be consis¬ 
tent with at least one pocket centered at F. For the 
purposes of determining the low-energy properties of the 
DOS, we therefore consider a single spherical pocket cen¬ 
tered at F. In this case, in the SCi state, there are two 
crossing nodal lines at = 0 and kp^x = 

SC2 phase, there remains a nodal line at kp^z = 0 and a 
couple of nodal points at kp^x = kp^y = 0 . Both nodal 
configurations are shown in Fig. [31 



Figure 3 . Nodes on a Fermi surface centered at F in the two 
different SC phases. Left: In the SCi phase, the Fermi surface 
contains the horizontal nodal line kz = 0 (red curve) and the 
vertical nodal line k^ + ky = 0 (blue curve). Right: in the SC2 
phase, the Fermi surface contains only one horizontal nodal 
line kz = 0 (red curve) and two nodal points (0,0, ikp) (black 
dots). 


We first consider the higher-temperature SCi state 
(T* < T < Tc). For concreteness, we assume Ai 0 
and A2 = 0, but the results are the same for the con¬ 


verse. The DOS is given by: 


p(w > 0 ) = ^ 6{uj - E{k)) 

k 



where Aq is the DOS at the Fermi level, and A((c) is the 
SC gap along the Fermi pocket centered at F. The latter 
can be conveniently described in terms of the polar and 
azimuthal angles 9 and (j) around the spherical pocket: 

A{k) = Ai sin — sm- ^—- 

Ri Ai sin20cos ^ , (10) 

where Ai has been rescaled by the factor Ai —>■ 
Ai -\/2 (kpa)'^ /8. There are two nodal lines at do = 7r/2 
{kz = 0) and (f)o = —f, ^ {kx = —ky ), which intersect 
along the Fermi surface at a single point, giving rise to a 
quadratic node. Expanding the DOS for small energies 
w ^ Ai, we find the asymptotic behavior: 


p(u;)«iVo(^log^ + |^log 4 ) (II) 

Fig. |T| compares the numerically-evaluated DOS with 
the asymptotic expression (EH) as function of energy, ev¬ 
idencing the dominant wlogw behavior at low energies. 


pi No pi No 



Figure 4 . The low energy DOS computed numerically (red 
curve) and via its asymptotic expression (blue curve). Left: 
DOS in the high-temperature SCi phase (T* < T < Tc). 
The blue curve is the asymptotic expression m with the 
dominant talogw behavior at low energies. Right: DOS in 
the low-temperature SC2 phase (T < T*). The blue curve is 
the asymptotic expression EH) with the dominant lo behavior 
at low energies. 

A careful analysis of Eq. (ITT]) reveals that while the 
linear-in-w term arises from the contribution of each indi¬ 
vidual nodal line, the w log oj term arises from the cross¬ 
ing point between the two nodal lines. To make this 
transparent, we expand the SC order parameter in the 
vicinity of one of the two crossing points, 9 o = 7 r/ 2 , 
^0 = 37r/4: 


A{9,cj)) = 2Ai{e-9o){^-M ■ 


( 12 ) 
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Clearly, the nodal dispersion is quadratic rather than 
linear— Substituting this expression in Eq. m and 
restricting the integration to the proximities of the cross¬ 
ing point, \9 - 9a\ ,\(j) - 4>o\ < A, yields: 


p{u}) Ki [ S [ d4>Re 

Jo Jo 


(13) 

where 9 = 9 — 9o, (j> = 4> — 4>o- A straightforward calcu¬ 
lation gives: 

p{uj) Ri log ■ (14) 

Clearly, the additional logarithmic contribution is a 
consequence of the quadratic dispersion near the crossing 
point of two nodal lines. Note that it will also give rise to 
logarithmic corrections to the low-temperature behavior 
of the thermodynamic quantities listed in the beginning 
of this section. More generally, tunneling experiments 
sensitive to p (ui) may in principle be able to identify 
this additional logarithmic contribution, which would be 
unambiguous evidence for this type of superconducting 
gap function. 

In the low-temperature SC 2 phase at T < T*, A 2 ^ 
0 and one of the nodal lines is replaced by two nodal 
points. As a result, the logarithmic corrections discussed 
above disappear. Indeed, the system in the SC 2 phase 
contains only one nodal line = 0 , and two nodal points 
kx = ky = 0. Computing the DOS at low energies, u <C 
A 2 , Ai, we hnd: 



p{uj) Ki No 


w , / 4Ai 

■ In ' 


2Ai y A 2 


U) 


4Ai A 2 


(15) 


As expected, the nodal line gives a linear contribution 
to the DOS, whereas the nodal points give a quadratic 
contribution. Of course, at energies higher than the 
scale of A 2 , but still much smaller than Ai, there is a 
crossover to the logarithmic behavior found in the SCi 
phase, which may be detectable experimentally. 

Interestingly, the coefficient of the linear-in-o; term de¬ 
pends on the ratio A 1 /A 2 , which in turn depends on the 
C 4 symmetry-breaking parameter ij, according to Eq. 0 . 
The dependence of this linear coefficient on p is shown in 
Eig.[5]below. This result reveals an interesting possibility 
to probe the impact of the tetragonal symmetry break¬ 
ing of the HO phase on the SC properties. Specifically, 
one could experimentally extract this linear coefficient - 
either from tunneling experiments or from penetration 
depth data - and study its dependence on compressive 
and tensile external strain, which would tend to increase 
or reduce the value of rj. 


pINg 



Figure 5. The linear coefficient of ui in the low energy DOS 
of the SC 2 phase as function of ^ — 1 oc »?. The red curve 
is the numerical result and the blue curve corresponds to the 
asymptotic expression m- 


III. ANISOTROPIC PROPERTIES OF THE 
SUPERCONDUCTING STATE 

Having established the thermodynamic properties of 
the SCi and SC 2 phases, we now discuss how the in-plane 
anisotropy appearing at the HO transition temperature 
is manifested in the superconducting properties. From 
symmetry considerations, one expects anisotropies in the 
in-plane penetration depth, critical magnetic field i?c 2 , 
etc. Experimentally, the behavior of the specific heat 
as function of the angle of an external magnetic field H 
has been used as a probe of the nodal structure of the 
superconducting gap in URu 2 Si 2 r^. Thus, given the ex¬ 
perimental feasibility and the existence of current data, 
we focus here on the magnetic-field angle-dependent spe¬ 
cific heat, C(^H), in both SCi and SC 2 phases. In the 
presence of an external magnetic field, a nodal supercon¬ 
ductor acquires a finite density of states p{H) at T — 0. 
As a result, the specific heat at low temperatures is given 
hy C{H)/T (X p[H). In a tetragonal system, p{H) dis¬ 
plays four-fold oscillations as the azimuthal angle of the 
applied magnetic field H changes. However, in an or¬ 
thorhombic system, which is the case of URu 2 Si 2 below 
the HO transition, two-fold oscillations are expected. 

A complete description of p{H), which is beyond the 
scope of this work, requires detailed knowledge of the 
Fermi surface and of the microscopic mechanisms of the 
superconducting state^i^. Here, we are interested in 
contrasting the anisotropies introduced by rj and man¬ 
ifested in p{H) in both SC phases, and also in study¬ 
ing the impact of twin domains present in real materi¬ 
als. Therefore, to keep the discussion general and mean¬ 
ingful, we consider a spherical Fermi surface and the 
semi-classical approximation for p{H^ first introduced by 
Volovik^. Valid in the regime of low fields and low tem¬ 
peratures, where the vortex lattice of the mixed state is 
dilute, this approximation focuses on the contributions 
of the extended quasi-particles around a single vortex, 
neglecting the contribution from the vortex core states. 
The main effect arises from the supercurrent Vs (r) gener- 
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ated by an isolated vortex, which causes a Doppler shift 
AE = k-Vs{^ in the quasi-particle excitation spectrum. 
Taking the spatial average over the unit cell of the vortex 
lattice then gives the zero-energy 


p[H) = 



(16) 


where A = $o /H = ttR^ is the area of the vortex lattice 
unit cell, which in turn is approximated by a circle of 
radius R = Here, $o is the flux quantum 

and + is the quasi-particle energy. For 

an arbitrary magnetic field along the h direction, H = 
Hfi, the supercurrent velocity Vg is given by: 


1 . , 

Vg =-n X r . 

2mr 


(17) 


To proceed, it is convenient to define the angle- 
dependent Doppler-shift energy: 


Eoik) 


kf 

2mR 


fix k 


(18) 


Evaluation of the DOS in Eq. (TIBI) then gives: 


P{H) 

No 


^ if EMk) < iA(k)i 

1 1 if En(k) > IA(k)l 


(19) 


It is convenient to define the ratio between the 
Doppler-shift energy scale and the magnitude of the SC 
gap: 


kp 

^ ^ 4mi7Ai ’ 


( 20 ) 


We consider first the SCi state, when Ai 7 ^ 0 and 
A 2 = 0. Since we are interested mainly on the in-plane 
anisotropies promoted by the HO tetragonal symmetry¬ 
breaking parameter rj, hereafter we consider the situation 
in which iJ is swept across the ab plane, and therefore 
characterized by the azimuthal angle (ph- Using the ex¬ 
pansion m for Ai yields the expression: 


This formula captures the general behavior of the 
numerically-calculated p{4>h), plotted in Eig. [51 In par¬ 
ticular, it reveals a clear in-plane two-fold anisotropy of 
p{4>h), with a maximum at (j)h = 7r/4, —37r/4 and a min¬ 
imum at 0ft = —7r/4, 37r/4. This behavior is a conse¬ 
quence of the quadratic node present in the Eermi sur¬ 
face at 00 = Tr/d, 9o = 7r/2 (see the previous Section): 
when E[ is parallel to the quadratic node momentum, 
the Doppler shift in the nodal quasi-particle spectrum is 
maximum, but when E[ is perpendicular to the quadratic 
node momentum, the Doppler shift is minimum. 





p(<P>,)lp(0) 



Figure 6 . The angle-dependent DOS in the presence of an in¬ 
plane magnetic field in the SCi state. The ratio between the 
Doppler shift energy and the gap is set to 7 = 0.1 in panel 
(a) and 7 = 0.01 in panel (b). 0 ft is the azimuthal angle 
of the magnetic field. The red curve corresponds to a single 
domain with A 2 = 0 , whereas the blue curve corresponds 
to a twin domain with equal-weight A 2 yf 0, Ai = 0 and 
Ai / 0, A 2 = 0. 


We now consider the SC 2 phase, which is the most rel¬ 
evant for comparison with experiments, since the latter 
are commonly performed at very low temperatures - pre¬ 
sumably below both Tc and T*. In the SC 2 phase, both 
Ai and A 2 are non-zero, but non-equal. We obtain: 


Eoik) _ Y^l - sin^ 6»cos2 (0 - 0ft) 


A{k) 


= 7 


|sin 20 cos (0 — f) I 


( 21 ) 


which can be substituted in Eq. (1191) for numerical eval¬ 
uation. In the weak field limit, 7 <C 1, we obtain the 
analytical expression: 


4Ao 


7 + 


COS 




( 22 ) 


Eoik) 


A{k) 


= 7. 


^1 — sin^ 9 cos^ (0 — 0ft) 


sin 29\ J cos^ ~ 7 ) + 




(23) 

which can then be substituted in Eq. (IT51) to compute the 
DOS. In Fig. [71 we plot p {(j)h) for two different values of 
the parameter 7 , considering the ratio A 1 /A 2 = 2. The 
behavior is similar to the SCi phase, namely, the DOS is 
two-fold anisotropic and maximum at 0ft = 7 r/ 4 , — 37 r /4 
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but minimum at (j)h = — 7 r/ 4 , 37 r/ 4 . The contrast be¬ 
tween the maximum and the minimum is proportional to 
the ratio A 1 /A 2 . In the Ci symmetric case, Ai = A 2 , 
to leading order in 7 , p{4>h) is angle-independent and 
constant, p{(l)h) = 47 / 7 r-|- 7 ^ ln( 2 / 7 ). 

An important issue is whether the two-fold anisotropy 
induced by the difference between Ai and A 2 (which in 
turn arises from the HO C 4 symmetry-breaking term t]) 
can be observed experimentally. In our analysis, so far 
we have considered only the ideal scenario in which a 
single C 2 domain is formed. In large samples - at least 
large enough to allow one to measure the specific heat 
- it is very likely that the system will break up in twin 
domains with rj > 0 (corresponding to A 2 = 0 in the 
SCi phase and Ai > A 2 in the SC 2 phase) and 77 < 0 
(corresponding to Ai = 0 in the SCi phase and Ai < A 2 
in the SC 2 phase). As a result, the two-fold anisotropy is 
washed out. To capture this effect, we calculated pt {(j^h) 
of a twin-domain system: 

Pt((('/i) = ^p((/'/i) + (</>/«+7r/2) • (24) 


P(tf'h)/P(0) 



p(<PhVm 



Figure 7 . The angle-dependent DOS in the presence of an 
in-plane magnetic field in the SC2 state (Ai,A2 ^ 0 ). The 
ratio between the Doppler shift energy and the gap is set 
to 7 = 0.1 in the upper panel and 7 = 0.01 in the lower 
panel. 4 >h is the azimuthal angle of the magnetic field. The 
red curve corresponds to a single domain with A1/A2 = 2, 
whereas the blue curve corresponds to a twin domain with 
equal-weight A1/A2 = 2 and A1/A2 = 1/2. The dashed 
black curve corresponds to the tetragonal symmetric system, 
with Ai = A2. 

The results for both SCi and SC 2 phases in the twin- 
domain case are shown by the blue curves in Figs. [5] 
and [T] While for the SCi phase there is a small four¬ 
fold anisotropy reminiscent of the anisotropies of the 


single-domain case, in the SC 2 phase the DOS is basi¬ 
cally angle-independent and indistinguishable from the 
tetragonal case (Ai = A 2 ). Interestingly, the in-plane 
angle-resolved specific heat data reported in^ is nearly 
flat as the field is swept from the a axis to the b axis. 


<n> 



Figure 8. Bulk averaged anisotropy (7) as function of the in¬ 
plain uniaxial strain e^y Without uniaxial strain, (7) = 0 
as averaged over different domains. It becomes nonzero upon 
application of uniaxial strain and remains non-zero even after 
the strain is released, inside the symmetry-broken phase. 

In order to decide whether the experimentally observed 
behavior is compatible with the orthorhombic SC 2 phase 
or with the tetragonal SC chiral phase, we propose to per¬ 
form angle-resolved specific heat experiments in samples 
under uniaxial external strain. In both cases, applying 
either compressive or tensile stress along the a (or b) axis 
will induce a two-fold anisotropy in p{4>h) ~ similarly to 
what we calculated for a single-domain case. This is not 
surprising, as the external strain explicitly breaks the C 4 
symmetry - in other words, it generates itself a term 7 in 
the free energy (| 3 ]) regardless of whether the symmetry 
was spontaneously broken at the HO transition. How¬ 
ever, upon releasing the strain, only in the SC 2 phase a 
non-zero 7 remains due to the alignment of the twin do¬ 
mains: This is nothing but the manifestation of the hys¬ 
teresis associated with the symmetry-broken phase (see 
Fig. [5]). Therefore, if the two-fold anisotropies of p{(j>h) 
persist after the applied strain is released, this would be 
unambiguous proof of the anisotropy in Ai and A 2 char¬ 
acteristic of the SC 2 phase. Note that devices capable 
of applying and releasing strain in a controlled way have 
been now widely used to study the nematic phase of the 
iron-based superconductors^. 


IV. COLLECTIVE MODES IN THE 
SUPERCONDUCTING PHASE 

As discussed above, inside the HO phase there are two 
SC transitions: at Tc, the system undergoes a normal- 
state to SC transition, characterized by the condensa¬ 
tion of one of the two components of the chiral SC or¬ 
der parameter (Ai if 7 > 0). At T* < a SC-to-SC 
transition takes place, in which the second component of 
the SC order parameter condenses (A 2 if 7 > 0) with a 
7r/2 relative phase with respect to the other component, 
breaking time-reversal symmetry. Detection of the sec- 





















ond SC transition in URu 2 Si 2 would be strong evidence 
for the scenario discussed here. Besides the usual ther¬ 
modynamic signatures of a SC transition, for instance 
in the specific heat, the system also displays distinctive 
collective SC excitations at T* - which in turn could 
be measured by spectroscopic techniques such as Raman 
scattering. 

To understand why this is the case, we note that 
the collective modes of a single-band superconductor 
are the phase and amplitude (also called Higgs) modes. 
While the former is always gapped by the Anderson- 
Higgs mechanism, the latter becomes soft near the su¬ 
perconducting transition. Near the usual normal-state 
to SC transition (such as the one taking place at Tc 
in our system), the energy gap of the amplitude mode 
is generally larger than the SC gap, implying that the 
mode decays into the particle-hole continuum. However, 
near a SC-to-SC transition (such as the one taking place 
at T* in our system) the situation is different, because 
the electronic spectrum has already been gapped by the 
hrst SC transition. As a result, it becomes in princi¬ 
ple possible to obtain a sharp SC amplitude mode at T* 
that does not fall into the particle-hole continuum^^— . 
The situation is similar to the appearance of the so-called 
Bardasis-Schrieffer mode in superconductors with nearly- 
degenerate superconducting states^^. 

To investigate this scenario, we include the time- 
dependence of the SC order parameter in the free energy 


A (A,) 


/i^+^|Ai|2 + ||A2|2 
f Af + f Af 
|A*A2 -f oAiA* 

V fA^A^ 


f A? + I Ai 
^^+u|Ai|2 + f|A2| 
|AiA2 

fAiA*+aAtA2 


Because we are interested only in the energy gap of the 
collective modes, and not in their dispersions, we ignore 
the spatial dependence of the SC gap. The energies of 
the four collective modes are therefore given by: 



(28) 


where \i are the eigenvalues of the matrix A. We first 
focus in the SCi phase, T* < T < T^, in which the 
equilibrium gap functions A^ g are given by Eqs. ([S)). Di¬ 
agonalizing the matrix A gives the following eigenvalues: 




Fsc 



dt 


= E7sc 
2=1 


aA, 

dt 


Fsc (Ai) 


(25) 


where Fsc is given by Eq. o and the coefficient of the 
time-dependent term is jsc > 0. Because we are in¬ 
terested in the behavior below F^, where the SC sys¬ 
tem is particle-hole symmetric, the time-dependent term 
must be quadratic in the time-derivative, in contrast to 
the usual linear iA*^ term that appears in the Gross- 
Pitaevskii equation, where particle-hole symmetry is not 
necessarily presen t^^^^ . 


To calculate the collective modes of the system, we 
linearize the gap equations near Tc by writing A^ = 
Ai_o + ^Ai, where Ai^ are the solutions of the time- 


independent equations, i.e. 



= 0. We obtain 


the coupled linear equations: 


/5A, 
I dAt 
dt'^ SA 2 
\6A* 


— A (Ai^o) 


/6A,\ 

SAl 

SA 2 

\SA*J 


with the matrix: 


(26) 


fAiA^+aA^Aa 

fA*A* 

!^+u\A2\^ + I\A,\^ 
f Af + f Af 


IA 1 A 2 \ 

|A*A2+aAiA* 
fAi + fA? 

^^ + w|A2p + ||Ai|V 


(27) 


Ai = 0 


A 2 

^3 


A 4 


—a + r] 

a {u + a — P) + r] {u — a + /3) 
2u 

a {u — a — P) + r] {u + a + P) 
2u 


with the corresponding eigenvectors: 


(29) 


vf = ( -1 1 0 0 ) 

= ( 1 1 0 0 ) 

= ( 0 0 1 -1 ) 

vj = ( 0 0 1 1 ) (30) 

Clearly, the two SC components Ai and A 2 decouple in 
this regime. The expressions for the eigenvectors reveal 
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that the eigenvalue Ai corresponds to the phase mode of 
the condensed SC component. Although in our calcu¬ 
lation it vanishes, it becomes a massive mode once the 
coupling to the electronic density is included (Anderson- 
Higgs mechanism)2ir— . The eigenvalue A2 is the ampli¬ 
tude mode of the Ai gap. Although it becomes soft at 
Tc , since Uc = rj, it does not give rise to a sharp collec¬ 
tive mode because at the electronic spectrum is barely 
gapped, implying that the mode falls into the particle- 
hole continuum. 

The other two eigenvalues correspond to modes of the 
incipient A2 component, which condenses only at T*. 
Because it condenses with a relative phase of 7r/2 with 
respect to Ai, A3 corresponds to an incipient amplitude 
mode, whereas A4 corresponds to an incipient relative 
phase mode. Because the latter does not couple directly 
to the electronic density, it does not become massive be¬ 
low T*—. Interestingly, we observe that at T*, which 
corresponds to: 


a* = —r] 


/u + P — a\ 
\u - /3 + a) 


(31) 


the mode A3 becomes soft. In this case, because Ai has 
already been condensed at T^, the electronic spectrum is 
gapped, implying that this collective excitation can be a 
sharp mode in the vicinity of T*. To analyze how these 
modes evolve below T*, we diagonalize the matrix A in 
the SC2 phase, T < T*, in which A^ 0 is given by Eq. 
0. We obtain: 


Ai = 0 


A 2 


—au 

u + fi — a 


A 3 


A 4 


—au 

u + /3 — a 
—2aa 
u + 13 — a 


I ^ I' u+ 13-a \ a^{l3-af- ' 

\u—j3 + a) {u + (3 — aY 

/ ^ + a^{j3-af ' 

y^ \u — l3 + a) {u + l3 — aY 

(32) 


Although the corresponding eigenvectors are straight¬ 
forward to obtain, we refrain from writing explicitly their 
lengthy expressions here. In Fig. we plot the eigenval¬ 
ues A2, A3, and A4 as function of temperature inside the 
SC state. As discussed above, A3, corresponding to the 
incipient A2 amplitude mode, vanishes at T*. This opens 
the interesting possibility of detecting the second SC 
transition spectroscopically. For instance, Raman scat¬ 
tering in the symmetry channel corresponding to the A2 
component (B^g/B2g irreducible representation of the or¬ 
thorhombic point group) could in principle detect a sharp 
mode near T* inside the gapped region of the spectrum 
(w < 2Ai). 

An important issue ignored in the analysis above is 
the fact that the superconducting state has nodal quasi¬ 
particle excitations, which in principle couple to the A3 


A 



Figure 9 . Eigenmodes of the SC state inside the HO phase, 
as explained in the main text: A2 (red), A3 (blue), and A4 
(green). Ai, corresponding to the global phase mode, is not 
shown, since it becomes massive due to the coupling to the 
electronic density. 


mode and can cause damping. To investigate this ef¬ 
fect, we compute the one-loop bosonic self-energy di¬ 
agram containing the coupling of A2 to the electronic 
states (and particularly to the nodal quasi-particles) at 
the vicinity of T*, where the A3 mode becomes soft (see 
Fig. (To]). At this temperature, the SC gap Ai is fully 
developed. Therefore, in Nambu space, the electronic 
Green’s function is given by: 


G{k,ivn) = 


dz A(k, 


27r ivn — z ’ 

with Vn = (2n -I- 1 ) ttT and the spectral function: 

= -^{Siz - Ej:) - S(z + Ej:)) 

X (zero - Ai(fc)cri) 


(33) 


(34) 


where Ag = ^^2-1-Af(fc) is the nodal quasi-particle 
excitation and ai are Pauli matrices in Nambu space. 
Thus, the bosonic self-energy shown in Fig. HOI becomes: 


n (icUn) 



dZ2 

27r 


^2 tr ^g(k)a2A (jt, zi^ g(k)a2A (jt, Z2J 


El 

k 


(iVm + iuJn - Zi) (Wrn “ Z2) 


(35) 


Here we set the external momentum q = 0 (since we are 
only interested in the dynamics) and included the Nambu 
vertex g(k)a2 = sin(20) sin((/) — 7r/4)CT2 corresponding to 
/S.2(k) = A2g{k)a2- A straightforward calculation gives: 


Il(iuJn) = 


d^k 

(27r)'^ 

1 


tanh 


f 3 E. 


ioj -b 2 E-^ iu! — 2 Ej: 


g\k) (36) 
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Figure 10 . Upper panel: Bosonic SC self-energy (in Nambu 
space) arising from the coupling between A2 and the quasi¬ 
particles near T*. Here g{k) = sin( 20 ) sin(<() — vr/l). Middle 
panel: the real and imaginary parts of the bosonic self energy 
when uj T. Lower panel: The ratio between the real and 
the imaginary parts of the frequency-dependent bosonic self¬ 
energy. It is clear that the real part is much smaller than 
the imaginary part at low frequencies. Moreover, since the 
imaginary part is proportional to u!^ and not to uj, the mode 
is under-damped. 


Performing the analytic continuation and subtracting 
the frequency-independent part, n(0), we obtain the 
imaginary and real parts of dll (w) = 11 (w) —11 ( 0 ) (here¬ 
after we consider w > 0 ): 


<5n"(a;) = TTtanh ( 


dn'(w) = 


d^fc 

(27r)'^ 


tanh 


• d^^fc 


S(uj - 2E^)g^(k) 

a^V(fc) 

El (4E| - 


(37) 


In the limit uj Ai, T, we find the low-energy asymp¬ 
totic behaviors of the bosonic self-energy: 


m"(w) 

(5n'(w) 


ttNq 

'W TAi 

No a;2 , f2A,\ 


(38) 


Thus, at low enough frequencies, not only does the 
imaginary part varies quadratically with the frequency, 
but also the real part is much larger than the imaginary 
part. Consequently, the A 3 mode can still be sharp near 
T* despite the damping introduced by its coupling to 
the nodal quasi-particles. To confirm these analytical 
results, in Fig. [101 we plot the behavior of dll" (w) and 
dn' (w) evaluated numerically from Eqs. (1571) . evidencing 
the sub-leading character of the imaginary part. 


V. CONCLUDING REMARKS 


In summary, we have investigated the impact of the 
tetragonal symmetry breaking promoted by the HO 
phase in the low-temperature chiral SC state of URu 2 Si 2 . 
Besides the anticipated splitting of the SC transition into 
two, the two resulting SC phases display very different 
low-energy behaviors. In particular, the nodal quasi¬ 
particle density of states of the higher-temperature SC 
phase acquires an anomalous logarithmic dependence due 
to the crossing of two nodal lines. Although absent in 
the lower-temperature SC phase at low energies, this log- 
behavior can in principle still be manifested for interme¬ 
diate energy ranges as a crossover effect. We have also 
shown the softening of one of the amplitude SC modes 
near the second SC transition, providing yet another 
signature of the interplay between tetragonal symmetry 
breaking and SC. Finally, we showed that the current 
angle-resolved specific-heat data is qualitatively consis¬ 
tent with either a tetragonal chiral state or an orthorhom¬ 
bic chiral state. We propose additional measurements in 
the presence of uniaxial strain to unambiguously distin¬ 
guish the two scenarios. 

It is important to critically analyze our results in face 
of recent data on the SC state of URu 2 Si 2 . First, our 
phenomenological model relies on the applicability of a 
Ginzburg-Landau approach. Although this seems to be 
the case in URu 2 Si 2 given the behavior of the thermo¬ 
dynamic quantities across Tc, large SC fluctuations have 
been recently proposed in Ref^S. As for the existence 
of two SC transitions, the temperature dependence of 
Hci has been interpreted as indirect evidence for one SC 
transition at Tc ~ 1.5 K followed by a second one at 
T* fv 1.2 K—. In contrast, recent Kerr data seem to be 
consistent with time-reversal symmetry being broken at 
T^. If this is indeed the case, it would imply that the 
tetragonal symmetry breaking at the HO transition is in¬ 
consequential for SC. On the other hand, the same data 
set reveals an anomalous Kerr signal well below Tc, at 
T Ri 1 K. This anomaly, combined with a “background” 
Kerr signal that onsets at high temperatures, indicates 
that at least for now one cannot rule out the possibility 
of two SC transitions in URu 2 Si 2 . As pointed out in Ref. 
[49], additional data are needed to settle this issue. An 
interesting possibility would be to perform Kerr measure¬ 
ments in strained samples. Due to the dependence of T* 
on the strain field, such a measurement would elucidate 

















11 


whether the anomalous Kerr signal could be a manifes¬ 
tation of a second SC transition. 

Our phenomenological results offer robust benchmarks 
that can be employed to study the interplay between the 
tetragonal symmetry breaking promoted by the HO state 
and the time-reversal symmetry-breaking promoted by 
the SC state in URu 2 Si 2 . With appropriate modifica¬ 


tions, our model should also be relevant to other systems 
in which chiral SC states have been proposed, such as the 
ruthenates^S and doped graphene^^. 
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